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Abstract 

The superstring and superbrane theories include gravity as a necessary and fundamen- 
tal part of a (future) unified field theory. Thus it is important to consider the alternative 
representations of general relativity as well as the alternative models of gravity. 
We study the coframe teleparallel theory of gravity with a most general quadratic La- 
grangian. The coframe field on a differentiable manifold is a basic dynamical variable. A 
metric tensor as well as a metric compatible connection is generated by a coframe in a 
unique manner. 

The Lagrangian is a general linear combination of Weitzenbock's quadratic invariants with 
free dimcnsionless parameters pi,p2,P3- Every independent term of the Lagrangian is a 
global 50(1, 3)-invariant 4-form. For a special choice of parameters which confirms with 
the local 5*0(1, 3) invariance this theory gives an alternative description of Einstcinian grav- 
ity - teleparallel equivalent of GR. 

The field equations of the theory is studied by a "diagonal" coframe ansatz which is a 
subclass of a most general spherical-symmetric Einstein-Mayer ansatz. The restricted La- 
grangian depends only on two free parameters pi, p^. 

We obtain a formula for scalar curvature of a pseudo-Riemannian manifold with a metric 
constructed from the static "diagonal" solution of the field equation. It is proved that the 
sign of the scalar curvature depends only on a relation between the parameters p\ and ps . 
Thus by a specific choice of free parameters a manifold of positive or negative curvature can 
be obtained. The scalar curvature vanishes only for a subclass of models with p\ = 0. This 
subclass includes the teleparallel equivalent of GR. 

We obtain the explicit form of all spherically symmetric static solutions of the "diagonal" 
type to the field equations for an arbitrary choice of free parameters. We prove that the 
unique asymptotic-flat solution with Newtonian limit is the Schwarzschild solution that 
holds for a subclass of teleparallel models with p\ = 0. Thus the Yang-Mills-type term of 
the general quadratic coframe Lagrangian should be rejected. 
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1 Introduction 



The resent progress in superstrings/M-branes theories and recognizing the supergravity as a 
low-energy limit of the supersymmetry renew an interest to alternative models of classical 
gravity. For instant, the generalization of AdS/CFT correspondence due to Witten and Yau || 
is based on manifolds with a positive scalar curvature. The metric of such manifold can not be 
obtained as a solution of the vacuum Einstein field equations which gives only the manifolds 
with zero scalar curvature. We will show in the sequel that the manifolds with a fixed sign of 
the scalar curvature are generated by generic solutions in teleparallel theories. 
The teleparallel description of gravity has been studied for a long time. The pioneer works of 
Cartan Q , Weitzenbock [Q] and Einstein [Q] dealt meaningful with the various models of unified 
(gravity-electromagnetic) field theory. 

The investigation in gauge field theory of gravity and in the Einstein-Cartan gravity (see ||, 
@],[|§ and the references) evoked considerable interest in teleparallel geometry. In the most 
general metric-affine theory of gravity (MAG) || the teleparallel Lagrangian appears as a 
separated part of the general MAG Lagrangian. 
For the recent investigations in this area see Refs. [|l0| to [p0| . 

The teleparallel approach to gravity can be viewed as a generalization of the classical Einstein 
theory. The main properties of this framework are as follows: 

1. The primary field variable of the theory (i.e. the coframe) has an intrinsic geometric 
sense. The basic geometrical quantities: metric, metric compatible connection and volume 
element can be constructed from the coframe in a unique manner. 

2. A covariant quadratic Lagrangian of the Yang-Mills type is defined in this teleparallel 
approach in contrast to classical general relativity (GR). 

3. A local energy-momentum current can be defined in a covariant manner. 

4. There exists a free-parametric wide class of possible field equations. So in order to get 
uniqueness additional symmetric requirements may be applied. For instance, the require- 
ment of a symmetry under the group of local ,50(1,3) pseudo-rotations uniquely yields 
the teleparallel equivalent of GR. 

5. The origin of gravity in this teleparallel approach is the energy- momentum current as well 
as the antisymmetric spinorial current. 

6. The more general modifications of the classical theory of gravity such as metric-affine 
gravity (MAG) include the teleparallel Lagrangian as a self-consistent sector. 

Let us give a brief account of teleparallel geometry. Consider an ra-dimensional differential 
manifold M endowed with a smooth coframe field {$ a , a = 0, ...,n — 1}. The 1-forms i? a is 
declared to be pseudo-orthonormal. This assumption completely determines a metric on the 
manifold M by the relation 

g = r,a b r®$\ (1.1) 
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where r\ ab = (+1, — 1, —1)- And conversely, the metric g is a unique metric that makes the 
coframe $ a pseudo-orthonormal. Thus it is possible to consider the coframe field i9 a as a basic 
dynamical variable and to treat the metric g as only a secondary structure. In order to have an 
isotropic structure (without peculiar directions) this coframe variable should be defined only 
up to global pseudo-rotations. In fact the true dynamical variable is the equivalence class of 
coframes with the pseudo-rotation as an equivalence relation. Thus in additional to invariance 
under diffeomorphic transformations of the manifold M the basic geometric structure has to 
be global SO(l,n — 1) invariant. Note that the coframe field is a complex of n 2 independent 
variables in every point of M while the symmetric metric tensor field has only n(n + l)/2 in- 
dependent components. 

An additional requirement of local SO(l,n — 1) invariance restricts the set of n 2 independent 
variables to a subset of n(n + 1)/2 variables, which can be related with the independent compo- 
nents of the metric. In this case the coframe (teleparallel) structure coincides with the metric 
structure. One can still describe this metric geometry by coframes, but in this case it does 
not represent a different geometrical structure but provides only a different mathematical tool 
- Cartan's "Repef e Mobile" . 

Thus teleparallel geometry can be considered as a certain generalization of the Riemannian or 
pseudo-Riemannian metrical geometry. In order to have a simplest (pure) teleparallel struc- 
ture one can choose the connections on the manifold M to be metric compatible. This metric 
compatibility of the connection is equivalent to the coframe compatibility: 

Dg = <==> Dtf a = 0, (1.2) 

where D is the covariant exterior derivative. 

Thus a 1-form of connection w a b is unique defined by the metric g and also, via the relation 



(1.1), by the coframe field "d a . Therefore the coframe field on the differential manifold completely 
defines a certain geometrical structure - coframe structure, which can be written symbolically 
as 

{M,r,g(r),w a b ^ a )}. (1.3) 

A straightforward generalization of this pure coframe structure can be obtained by considering 
the connection w a b as an additional independent primary variable. Such generalization, which 
we write symbolically as 

{M,r,g(r),w a b }. (1.4) 
can be refereed to as Cartan teleparallel geometry. 

The structure can further be generalized by taking the metric field g as an independent geo- 
metrical variable: 

{M,r,g,w a b }. (1.5) 

This most general geometrical structure is realized in the framework of metric-affine gravity 
(MAG) §• 

In the present work we study the general teleparallel model of the type ( |1.3[) with a quadratic 
Lagrangian on a 4D-manifold of Lorentzian signature {1,— 1,— 1,— 1}. 



3 



2 Quadratic Lagrangian 



The general quadratic (in the first derivatives of the coframe field $ a ) Lagrangian can be written 
as a linear combination of the Weitzenbock second order invariants Q. This Lagrangian is 
known as a translation invariant Lagrangian of Rumpf E| : 



1 3 

where I is the Plank length and * is the Hodge dual operator. [] 

The Lagrangian includes the strength d$ a and its irreducible (under the global Lorentz group) 
parts ^d'&a. The explicit expressions for these irreducible parts are: 

Wdtf" = d$ a - (2 U$ a - ( 3 W a , (2.2) 

( 2 W a = -i9 a A (e;,Jdi? 6 ), (2.3) 
3 

m M a = i e «j(0& A( J0 6 ). (2.4) 
o 



The coefficients 01,02,03 in ( |2.l[) are dimensionless constants. Rearranging the different parts 
of the Lagrangian (^]l]) results in a linear combination of quadratic terms 

i=l 

with 

(1 V = dtf a A*d$ a , (2.6) 

( 2 ) y = (d4Ar)A*(4A#), (2.7) 

( 3 ) y = ^ a A/)A*(^Ai?j. (2.8) 
The dimensionless coefficients en and pi are related |18| as: 



111 

Pi = g (02 + 2ai) , P2 =g (03-01), P3 = g(ai-a 2 ), (2.9) 

01 = pi + P3 , o 2 = pi - 2p 3 , a 3 = pi + 3p 2 + p 3 • (2.10) 



Every independent part of the Lagrangians ( |2,lD and (|2.5| ) is diffeomorphic covariant and in- 
variant under the global (rigid) Lorentz group - translation invariant. Thus different choices of 
the free parameters or pi yield different translation invariant models of gravity. 
Let us list the o, and the pi coefficients for different teleparallel Lagrangians [18]: 

1 Here and later the down-indexed coframe is obtained by the Lorentzian metric, thus i?° = #o> = — 
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Note also the following possibility to modify the Lagrangians (2.1) and Q2.5p . One can add to 
them a zeroth part 



n = -M a A *0 a . 
4 



(2.11) 



It is easy to see that this term is an analog of the cosmological A-term in Einsteinian gravity. 
We will not use this term in the sequel. 



3 Field equation and ansatz 

In order to obtain the field equation for the coframe field t? a one has to calculate the variation 
of the Lagrangians ( |2.1| ) or ( |2.5| ) relative to a free variation of the coframe. This procedure is 
not quite regular because of the Hodge operator *, which itself depends on the coframe field # a 
and should also be a subject of variation. Thus we have, in general, for an arbitrary form w, a 
non-commutative relation 

5 * w — *5w / 0. (3.1) 
The variation of the Lagrangians (2.1), (2.5) can be produced in the following ways: 



1) Directly by using the Euler-Lagrange equations |2 

/ dV \ dV 

2) By using the master formula of [18] - the explicit expression for the commutator ( |3.1| ). 

3) By using a formula [21] for variation of the quadratic Lagrangians of the form 5(a A */?). 
All these different methods of variation result in the same second order field equation, that we 
will write symbolically as 

px (2«L a + 0) Qa _ 2 Wg tt ) + P2 ( - 2^L a + 2( 3 ) Q a + ^Q a - 2^Q a ) + 

P3 ( - 2^L a + 2^Q a + ( ? )Q a - 2^Q a ) = 0, (3.3) 

where the leading part (which includes all the pieces with the second order derivatives) is a 
linear combination of the terms: 



«L a = d*dtf OJ (3.4) 
^L a = tf a Ad*(cft? 6 Atf 6 ), (3.5) 
(3) L a = A(i* (tf a AcW b ). (3.6) 
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The quadratic part (which involves only the pieces quadratic in the first derivatives) is a com- 
bination of the terms 



(1) Qa 


= e a _\(d$ b A*d$ b ), 


(3.7) 


{2) Qa 


= {e a AM b ) ^*d^ b l 


(3.8) 




= dtf a A*(dtf b Ai?b), 


(3.9) 


(4) Qa 


= e a J (cM c A i? c A *((M h A 6 )) , 


(3.10) 


(5) Qa 


= (e a JA? fe ) Atf b A*((M c Atf c ), 


(3.11) 




= A*(tf a A(M fe ), 


(3.12) 


(7) Qa 


= e a J (tf c A cM 6 A *(d$° A l? 6 )) , 


(3.13) 


(8) Qa 


= (eaJdi? 6 ) At? c A*(cM c A0&). 


(3.14) 



We are mostly interested in static, spherical symmetric solutions of the equation (| 
well known that a general static spherical symmetric metric can be written as 



ds z 



^dt 2 



,2.g 



(dx 2 + dy 2 + dz 2 ), 



where / = f(r), g = g(r). 

Due to Einstein-Mayer H a general static, spherical symmetric coframe has the form 



if' 



-Jdt, 



if = wx { dt + eW, 



1,2,3, 



where / = f(r), g = g(r) and w = w(r). 

We begin the consideration with a simplified "diagonal" ansatz 

^0 = e f(x,y,z) d ^ = e g(x,y,z) dx i^ • 



1,2,3. 



It is 



(3.15) 



(3.16) 



(3.17) 



and in a sequel we will study the spherical-symmetric "diagonal" solutions to the field equation 
( |3.3| ). Thus we restrict our attention to a sub- family of solutions with an identically vanishing 
function w. This restriction is mainly made in order to simplify the calculations. The con- 
sideration of a non-restricted spherical-symmetric coframe as well as the study of the physical 
sense of its "non-diagonal" part will be made in a sequel paper. Note, however, that even the 
restricted coframe leads to a most general spherical-symmetric metric (3.15). Thus it can be 
believed that the "non-diagonal" pieces of the coframe can be related to an additional (possible, 
non-gravity) interaction. 

For the "diagonal" coframe ( |3.17| ) the following "parallel" relation holds 



a A d$ a = 0. 



(3.18) 



This condition is diffeomorphic and global SO (1,3) invariant, thus it can serve as a necessary 
condition to have a "diagonal" coframe. The relation ( |3.18 ) shows that the multiplier of the 
coefficient p2 in the Lagrangian ( |2.5) ) is identically zero. Therefore the field equation fl3.3| ) 
reduces to: 



Pi 



2 (1) L a + «Q a - 2^Q a ) +pJ- 2®L a + 2^Q a + ^Q a - 2^Q a ) = 0. (3.19) 
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The coframe ( [3.17D being substituted (Appendix A) in the field equation ( 3.19| ) splits it in the 
temporal (a = 0) and in the spatial (a = i = 1, 2, 3) parts. The temporal part takes the form 



Pl ( - 2A/ - 2(V/ • V 5 ) + 2{Vg) 2 - (V/) 2 J * + p 3 ( - 4A 5 - 2(V<?) 2 J * = 0, (3.20) 
which is equivalent to a scalar equation 

P! (2A/ + 2(V/ • V<?) - 2(Vg) 2 + (V/) 2 ) + 2p 3 (2Ag + (V 5 ) 2 ) = 0. (3.21) 
As for the spatial part of the field equation ( |3.19| ) 

px M 2A 5 + 2(V/ • V 5 ) - (Vf) 2 )rjik + 2g ik + 2f m - 2fif k - 2g i9l \ * i?* + 

2p 3 ((A/ + A<? + (V/) 2 )^ + /« + + fif k - fm - 9i f k - g igt ^ *ti k = 0. (3.22) 

The system ( |3,20D and Q3.22| ) is exhibited in the scalar form by 

Pl f 2A/ + 2(V/ • Vg) - 2(Vg) 2 + (V/) 2 ) + 2p 3 (2Ag + (Vg) 2 ) = 



pi ( (2A<7 + 2(V/ • Vg) - (V/) 2 ) r] lk + 2 5ifc + 2f m - 2f t f k - 2g i9} \ ( 3 . 23 ) 

+2p 3 H A/ + Ag + (V/) 2 ^77ife + /ifc + ftfc + /j/fe - - ffi/fc - ftfite^ = 0. 

The last equation is symmetric in respect to transposition of the indices (i <—> j). Thus the 
result is an over-determined system of 11 equations for two functions / and g. This system 
certainly has a trivial solution / = g = which yields a flat Minkowskian metric. The interesting 
fact is that this over-determined system has also nontrivial solutions for an almost all values of 
the free parameters pi,P3- The explicit form of these solutions will be obtained in the sequel. 



4 Curvature 

In this section we establish a formula for scalar curvature of a manifold, which metric is con- 
structed from the solution of the system (3.23). The result does not depend on a specific 



solution of the field equation. For actual calculation of the scalar curvature we use the formulas 
p9| for non-vanishing components of Ricci tensor in the case of "diagonal" metrics. 
Let the components of the metric tensor be 



g aa = e a e 2Fa , g af3 = for a ^ /3, (4.1) 

where e Q = ±1. 

The components of Ricci tensor are: 

R a p = Yl { F hP F P,<* + F <*. A* ~ F„, a F„,p - F^aj) for a ^ 0, (4.2) 

tM^a,P 



J In all formulas in this section the summation over repeated indices is not used. 
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TP TP TP Z TP 

L a,ot-L fi,a 11,01 A^QV 



- E[ 

+ P P (ZiFcc-F^) ( p jp _ jp2 _ jp _ Tp ST^ vp 



For the diagonal coframe ( 3.17] ) we have 

e = -1, ei = 1, F = f, Fi = g for i = 1, 2, 3. 
Hence using (|4.3| ) we obtain 



. e 2(/- 9 ) ( (F^Fo,! - Fq 2 ! - F ,i,i - F 0i1 (F 2j1 + F 3il )j + 

^2,2-^0,2 - Fl 2 - F , 2)2 - F , 2 (iq,2 + ^3,2)) + 



-^3,3-^0,3 — Fq, 3 ~~ -^0,3,3 _ -fo,3(-^l,3 + -^2,3, 



= 2(/- 5 ) { f 2 



ft + hi + fl9l + ft + Hi + /252 + ft + /33 + /353 

Thus the temporal component of the Ricci tensor is 

R 00 = e W-9) ( A / + (V/) 2 + (V/ • Vg) 

Using (4.2) we compute 



R 



ii : 



^1,1-^0,1 — ^0,1 ~~ ^b,i,i + eieoe 2 ^ 1 ^Fo.o-^i.o — ^f,o ~~ -^1,0,0 — -^1,0(^2,0 + -^3,0, 
^1,1-^2,1 — -^2,1 ~~ -^2,1,1 + eie2e 2( - Fl F ^ ^2,2-^1,2 — F\,2 ~~ -^1,2,2 — -^1,2(^0,2 + -^3,2, 
^1,1-^3,1 — -Pf,i — -^3,1,1 + eiese 2 ^ 1 F ^ (V33.F13 — F 2 3 — ^1,3,3 — i 7 !, 3(^0,3 + -^2,3, 

5ii + (533 + 9z{h +93) 



9i h ~ fi ~ hi ~ 9ii + [922 + 52(/2 + 92) 
Thus 

Ru = -hi - 911 + 2gih - ft + g\ -&g- (Vg • V/) - (Vg) 2 . 

Therefore, the spatial components of the Ricci curvature tensor are 

Ru = -hi - 9» + 2gih -ff + 9t - As - {Vg ■ Vf) - (Vg) 2 . 
Consequently, the scalar curvature takes the form 



R 



e -2 9 Uf + 2Ag + ( V/ )2 + (V/ . Vg) + {Vg) 2\ _ 



(4.3) 



(4.4) 



(4.5) 



+ 



(4.6) 



(4.7) 



(4.8) 



Note that this formula is valid for the diagonal coframe ( 3.17 ) even if it does not satisfy the 



field equation (|3.23| ). Calculate now the scalar curvature of the pseudo-Riemannian manifold 



constructed from a solution of the field equation. Taking the trace of the second equation of 
the system ( |3.23| ) and adding the first equation we obtain 

2 P1 (2A 5 + A/ + 3(V/ • Vg)) + 4p 3 (2Ag + A/ + (V/) 2 + (V/ • V<?) + (V 5 ) 2 ) = 0. (4.9) 

Substituting the second derivatives expression from the formula ( |4.8| ) we get 

Re 2 9( Pl + 2p 3 ) = - Pl (Vf - Vg) 2 . (4.10) 

We will see later that in the special case p\ + 2/?3 = the system ( |3.23| ) has a non-trivial solution 
only if / = g. 

Thus the sign of the scalar curvature depends only on the relation between free parameters 
pi,P3 and does not changes on the manifold. This sign is an invariant for every particular 
teleparallel model (particular choice of the parameters pi,p 3 ). 
The conclusion is: 

Proposition 4.1: Suppose that p\ + 2p 3 ^ 0. The "diagonal" solutions of the field equation 



(3.S) describe the pseudo-Riemannian manifold with the scalar curvature 

(4.11) 



R= £L^ e - 2 °(Vf-Vgf 

Pi + 2/9 3 



The sign of the scalar curvature R is opposite to the sign of p\{pi + 2p 3 ) and it is an invariant 
of a particular teleparallel model. 

The scalar curvature is zero if and only if p\ = 0. - the teleparallel equivalent of GR. 

The case p\ + 2p 3 = will be treated later. Note that the subclass p\ = includes the Einstein 
theory of gravity in the form of teleparallel equivalent of GR. 

We have shown that the scalar curvature of manifold constructed from a "diagonal" solution of 
a specific teleparallel model has a fixed sign. In order to clarify the result let us recall a similar 
fact from the geometry of surfaces. Consider, as an example, a 2D-surface with a metric 

dl 2 = e 2f {du 2 + dv 2 ), (4.12) 

where / = f(u,v). Note that for an arbitrary 2D-manifold every metric can be transformed to 
this conformal invariant form. The Gauss curvature for the metric ( |4.12| ) takes the form 

K = -e' 2f Af. (4.13) 

Consider now a class of metrics that satisfy a "covariant field equation" : 

A/ = p(V/) 2 . (4.14) 

The Gauss curvature takes now the form 

K = -p(e- 2 f{Vf) 2 ). (4.15) 



The term in the brackets is positive, thus the sign of the curvature depends only of a free 
parameter p of the model. Note that fl4,15| ) is closely similar to the formula ( 4.11 ). The 
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curvature vanishes only in the case p = 0, i.e. / is a harmonic function. This is a 2D Euclidean 
analog of Einstein gravity. 

In fact we have proved this property of having fixed sign only in a case of "diagonal" metric. 
Our conjecture that this property can be proved for an arbitrary static metric that satisfies the 
general teleparallel field equation. 



5 Spherical-symmetric ansatz 

Consider a solution of the field equation which depends only on one radial variable s = x 2 + 
y 2 + z 2 . 

&° = e K s Ut , ^ = e 9(s) dt\ i = 1,2,3. (5.1) 

With this ansatz we reduce the general problem of determining 16 unknown functions (the 
components of the coframe i? a ) of four variables (coordinates) to determining two functions / 
and g of one radial variable s. We obtain, instead of ( |3.23[ ), an over-determined system of three 
second order ODE for two independent variables 

' p x (2f"s + 3/' + 2f'g's - 2(g') 2 s + (f) 2 s) + 2p 3 fa's + 3g> + (g') 2 s) = 
< Pl (2g» + 2/7 - 2(f) 2 - 2(g') 2 ) + 2p, (/" + g" + (f) 2 - 2fg> - (g') 2 ) = (5.2) 
_ p x (4g"s + 4g> + Afg's - 2(f) 2 s) + 2p 3 (2f"s + 2f + 2g" s + 2g' + 2(f) 2 s) = 0. 

We will see later that for certain special values of the free parameters pi, P s these three equa- 
tions are algebraically dependent and the system fl5.2j ) reduces to a good posed system of two 
independent equations for two independent variables. For most of the values of free parameters 
Pi,P3 the three equations of (|5.2| ) are independent. Using this fact we obtain an explicit form 
of exact solutions. Let us make an algebraic transformation of the system ( ^.2|), which allow to 
decrease the order of the system. The subtraction of the second equation of (fT^) multiplied by 
2s from the third equation multiplied by 3 results in 

Pi ( V'a + 6g' + 4/7 « - (f) 2 s + 2(g') 2 s) + 

2 P3 (2f's + 3f + 2g"s + 3g' + 2(f) 2 s + 2/7* + (</) 2 *) = 0. 

The subtraction of the second equation of ( |5.2p multiplied by (2s) from the third equation 
results in 

Pl + (f) 2 s + 2(g') 2 s) + 2^ (/' + g' + 2fg's + (g') 2 s) = 0. 
Thus we obtain an equivalent system to ( |5.2j ) 

' Pl {2f's + 3/' + 2/7* - 2(g') 2 s + (f) 2 s) + 2p 3 fa's + 3g' + (g') 2 s) = 
p 1 (4g"s + 6g' + Afg's - (f f s + 2(g') 2 s) + 

2p 3 (2f's + 3f + 2g"s + 3g' + 2(f) 2 s + 2fg's + (g') 2 s) = 5 ' 3 
^ Pl (2g> + (ffs + 2(g') 2 s) + 2p 3 (/' + g> + 2fg's + (g') 2 s) = 



10 



We can now decrease the order of the system (^^) by introducing new variables 

u = j , v = g S2 . 
WE obtain a system of first order ODE 



(5.4) 



pi \2s2u' + 2uv - 2v 2 + u 2 j + 2p 3 \2s^v' + v 2 j = 

pi Us 2V + Auv - u 2 + 2u 2 ) + 2p 3 ^2slu' + 2slf/ + 2u 2 + 2uu + v 2 ^ =0 

pi ( 2us5 + u 2 + 2u 2 ) + 2p 3 (us^ + us2" + 2uu + v 2 ) = 0. 



(5.5) 



Note that the first two equations of ([Si]) contain the first order derivatives and the quadratic 
expressions of the "strength" variables u and v. This is a general form of a non- linear (Einstein 
type) equation in gravity. The third equation of ([5j]) is an algebraic relation between the 
variables u and v and the radial coordinate s. 

Solving the first two equations of ( |5.5[) for the derivatives terms results in 



(2p 3 + Pi)(ps - Pi) [2s2u'j=(p 2 1 + 2p 3 pi - \p\)u 2 - 2{2p\ + pip 3 - p 2 )uv - 2pi(p 3 + Pl )v 2 
2(2p 3 + pi)(p 3 - p 1 )(2sh') =pi(2p 3 - pi)u 2 + 4p 2 uv - 2(2p\ - 3pip 3 - p\ )v 2 



pi I 2vs 2 + u + 2u ) + 2p 3 1 us 2 + vs2 + 2m; + v 



0. 



(5.6) 



The three equations (5.5) or ( |5.6j ) contain the derivatives u',v' so one can hope to obtain an 
algebraic equations for the functions u and v. In order to obtain such algebraic relation we 
writhe the second and the third equations of the system ( |5.5| ) as 

(5.7) 
(5.8) 



(piv + psu + r 3 v') = Fs 3 / 2 , 
(piv + p-3U + r 3 v) = Gs~ 1/2 , 



where the quadratic functions of u and v are 



F 
G 



F(u,v) 
G{u,v) 



-/'))«"+ I (/'I +/>:■!)'"• + 2 (pj ~P3)r 

1 



- pin + 4p 3 ui; + 2(p - 1 + p 3 )?/ 



Taking the derivative of the ( |5.8] ) and comparing it with (5.7) we obtain 

dG . dG . ( „ 1 _\ i 

The derivative terms in ( |5.6[) can be written symbolically as 

u' = (j>i(u, v)s~ 3 / 2 , v' = 2 (u, v)s~ 3 / 2 . 



(5.9) 
(5.10) 

(5.11) 
(5.12) 



Substituting this equations in ( |5.11| ) we obtain an algebraic equation that steel contains the 
independent variable s 



dG ± dG ± / 1\ 1/2 
-5-^1 + ~^-<h = (F+-G)s 1/2 



du 



dv 



(5.13) 
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Deriving now s from fl5.8| ) and substituting in ( f>. 13 ) we obtain an algebraic equation for u and 
v which we are looking for 

d^ 1+ d^ 2 ) = ( F+ 2 G ) G ' (5 ' 14) 

In order to write this equation in the explicit form we make a use of an algebraic computations 
program "REDUCE" . The result is 

Pi(pi - Pz){u ~ v)(piv + p 3 u + p 3 v) (piu 2 + 4p 3 uv + 2 (pi + pz)v 2 ^ = 0. (5.15) 

This equation gives us a necessary condition to obtain a solution of the system (|5,5| ) or ( |5.6|) . 

Proposition 5.1: In order to have a nontrivial spherical- symmetric static "diagonal" solution 
of the field equation (3.5) one of the following condition should be satisfied: 

1. The free parameters satisfy the relation 

Pl(pi-p 3 ) = 0. (5.16) 



2. Function u is zero thus 



3. Function v is zero thus 



Pi + P3 = 0. (5.17) 



P3 = 0. (5.18) 

Jf.. The nonzero functions u and v satisfy the proportional relation 

u = Xv, A = A(pi, p 3 ). (5.19) 



6 Solutions 

Now we are able to derive the explicit form of the solutions for the system ( |5.2|) for different 
values of the free parameters pi,pz- It is useful to write the formula ( 4.11| ) for the scalar 
curvature via the radial functions u and v. Using the relations 

Vf 2 = A(f'fs, Vg 2 = 4(g') 2 s, and V f ■ Vg = Af g' s. 

we obtain 

R = ^-e-^(Vf-Vg) 2 = -4^—e- 2 ^(f-g') 2 . 

Pi + 2p 3 pi + 2p 3 



Thus the scalar curvature is 



R = -- 2 ^-e- 2 *(u-v) 2 . (6.1) 
s z px + 2p 3 
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6.1 GR-type Lagrangian: p\ = 



The system ( |5.5[ ) takes the form 

2s2 v + v 



2s§*/ + 2s §t/ + 2-u 2 + 2uv + v 2 = 



(6.2) 



lis 2 + vs 2 + 2-u-u + v = 
Integrate the first equation to obtain 



Cy/s + V 



(6.3) 



where c is a constant of integration. 

Substituting this relation in the third equation of the system (|6.2[), we find 



cs 



c 2 s — 1 



(6.4) 



Substituting the functions (6.2 - 16.4| ) in the second equation of fl6.2|) we obtain an identity. 
Integrate Fl the functions u and w to obtain the ansatz functions (r = y/s): 



f = ln 



1 



CT 



1 + 



1 ' 



5 = 21n(l + -). 

cr 



(6.5) 



By taking the parameter of integration to be inversely proportional to the mass of the central 
body c = — we obtain the coframe field in the form 

•J m. 



tf 



i m 

1 l TJL > 

1 2r 



1?' 



(1 + ^) 2 ^, 



1,2,3. 



This coframe field yields the Schwarzschild metric in isotropic coordinates 



ds 1 



1 



m\ 2 
2r 



(it 2 



l + ^) 4 ((ix 2 + (iy 2 + dz 2 ). 



(6.6) 



(6.7) 



1 + — 

Note that the parameter p2 is not determined via the "diagonal" ansatz. Thus the Schwarzschild 
metric is a solution for a family of teleparallel field equations which defined by the parameters: 
pi = and p2 , P3 - arbitrary. Note that for the teleparallel equivalent of GR we need beside of 
zero p\ an additional condition: p^ + 2p2 = . 



6.2 pi = p 3 

The system ( |5.5[ ) takes the form 

2s§ (v! + 2v') + u(u + 2v) = 

4sf (u' + 2v') + (n + 2t>)(3n + 2t>) = (6.8) 
2s^{u + 2v) + (u + 2t>) 2 = 0. 

3 The constant of integration on this step can be omitted, because we can still re-scale the coordinates. 
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Subtracting the first equation (multiplied by 2) from the second one we obtain that the system 
is satisfied if and only if 

u + 2v = 0, ==> f + 2g = 0. (6.9) 

Consequently the coframe is 

^0 = e ~2g(r) d ^ 0i = e 9(r) dx i ( g 1Q) 

and the metric is 

ds 2 = e~ i9ir) dt 2 - e 2 ^ r \dx 2 + dy 2 + dz 2 ). (6.11) 

Thus the solution incorporates an unknown function g = g(r). The scalar curvature of the 
metric ( 6.11| ) is negative for an arbitrary choice of the function g = g(r). 

R = -3(g') 2 e- 2g . (6.12) 



6.3 u = ==> px + p 3 = 

The equation ( |5.15 ) can have a non-trivial solution with a zero function u if and only if the 
parameters satisfy 

Pi + P3 = 0. (6.13) 

This is a case of FAff La grangian 

C = d^ a A *d^ a , (6.14) 

where d^ is a coderivative operator. 

The second and the third equations of ( |5.5[ ) are satisfied identically for zero u and parameters 
satisfied (|6.13| ). As for the first equation of (5.5) it gives for function v 

sl v ' + v 2 = 0. (6.15) 

The solution of this equation is 

-■ ^ => 9 ' 1 



Cy/s-2 y s{cyfs-2)' 

Omitting the constant of integration we obtain the solution 

/ = 0, g = In (l - -) = In (l + ^ 

\ crJ V r 

where we introduce a new length dimensional parameter tq. 
Thus the correspondent coframe is 



= & = (1 + — )dx\ i = 1,2,3. (6.16) 



r 
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This coframe generates an asymptotically-flat metric 

ds 2 = dt 2 - (l + ~Y (dx 2 + dy 2 + dz 2 ). (6.17) 

Let us first determine the sign of the parameter tq. The metric ( |6.17| ) represents a point-like 
solution if the correspondent ADM- mass is positive. Rewrite the metric ( |6.17 ) in the spherical 
Schwarzschild-type coordinates (where the circumference of a circle with center at the origin is 
equal to 2nr). In the spherical coordinates the metric is 

ds 2 = dt 2 -(l + r -±)\dr 2 + r 2 dtf). 

Using now the translation 

r = r + tq 

we obtain the asymptotic-flat metric in the Schwarzschild coordinates 

ds 2 = dt 2 - df l . 2 -r 2 dn 2 . (6.18) 
The ADM mass for the metric ( |6.18| ) takes the form 

m := Jim T - (\ - (1 - y ) 2 ) = r . (6.19) 



r— »oo 



Thus by taking the parameter ro to be positive we obtain a particle-type solution with a finite 
positive ADM-mass. 

The metric ( |6.17 ) is singular at the origin of coordinates r = and consequently the metric 
( 6, 18] ) is singular at r = ro. In order to clarify the nature of this singularity compute the scalar 
curvature of the metric ( 6.17| ) via the formula ( |6.1| ). The result is 



r 2 
' n 



R = r J im = ^- (6 - 20) 

(ro + r) 4 r 4 

This function is non-zero and regular for all values of r including the origin (r = 0). 

The proper distance for a radial null geodesic in the metric ( 6.1S| ) is equal to the proper time 

and attach the infinity 

T 1 df 
l = t= / TTa^ 00 

Thus the point r = (f = ro) does not belong to any final part of the space-time. 



6.4 v = ==> p 3 = 

The equation ( |5.15| ) has in this case a non-trivial solution only if p% = 0. The system ( |5.5D 
yields u = thus it does not has a nontrivial solution. 
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6.5 Proportional solution. 



We see from the equation ( |5.15 ) that generic solutions u and v satisfy the homogeneous algebraic 
equation thus they should be proportional one to the other 



u = Xv, 

Rewrite the third equation of (|5.5|) as 



A = A(pi,p 3 ). 



U V 



Substituting u = Xv we obtain an algebraic linear equation for (^-^j with constant coefficients. 
Thus the general solution of the system should be of the form 

U = dy/s, V = by/s, 



where a and b are constants that depend on the parameters pi,p 3 . 
Consequently the ansatz functions / and g satisfy Q 



r 



9 



us 



vs 



-3/2 _ « 

s 

-3/2 = b 



> 



> 



f = aln-y, 



g = bln-y. 



The scalar curvature (|6.l| ) for this type of solutions takes the form 

, pi rt b {a - bf 



R 



Pi + 2^3 



s 16+2 



(6.21) 

(6.22) 
(6.23) 

(6.24) 



Substituting the solutions ( |6.21 ) in the system ( |5.5| ) we obtain an algebraic system for two 
scalars a and b: 



pi (a + 2ab - 2b 2 + a 2 ) + 2p 3 (6 + b 2 ) = 

pi (26 + 4a6 -a 2 + 2b 2 ) + 2p 3 (a + b + 2a 2 + 2a6 + 6 2 ) = 

pi(26 + a 2 + 26 2 ) + 2p 3 (a + b + 2ab + b 2 ) = 0. 



(6.25) 



Extracting the second equation from the first and from the third we obtain an equivalent system 

' pi(a + 2ab - 2b 2 + a 2 ) + 2p 3 (6 + 6 2 ) = 
(a + 6)(pi(26- a) + 2p 3 a) = ( 6 .26) 
a(pi(26- a) + 2p 3 a) = 0. 



The case a = and consequently u = was treated above thus we obtain from the third 
equation of this system that the nontrivial solutions should satisfy 

pi(26 -a) + 2p 3 a = 0. 



4 By rescaling the coordinates in the coframe and in the line element we can use the same constant of integration 
in the functions / and g. 
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Rejecting the case p\ = (also treated above) we obtain 

Pi ~ 2/93 



b = a- 



2pi 



Substituting this relation in the first equation of ( |6,26| ) we obtain 

(Pi - P3)(3pi - 2p 3 )(pi + 2p 3 )a = 2pi(p 3 - pi)(pi + 2p 3 ). 
In order to satisfy this equation we have to separate the following three cases: 

• Pi = P3- 

Via ( 6,27 ) the constants are connected as 

a = -2b 



and the system (6.26) satisfied identically. 
Hence the coframe 



if 



?•() 



2b 



dt, ffi 



includes two free parameters ro and b. The correspondent line element is 



ib 



dt 2 



2b 



(dx 2 + dy 2 + dz 2 



pi = -2p 3 



Via ( 6.27 ) the constants in this case are equal 

a = b 



and the system (6.26) satisfied identically. Hence the coframe is 



= ( — ) dt, & = ( — ) dx\ 



T \ a 



f \ a 



ro 



The correspondent line element is conformal flat 

. 2a 



ds 2 



(—) 2a (dt 2 - dx 2 - dy 2 - dz 2 ). 
Vrn/ 



(6.27) 



(6.28) 



(6.29) 



(6.30) 



(6.31) 



(6.32) 



• 3pi = 2p 3 . 

The equation ( |6,2S{ ) does not satisfied for any finite value of a thus the field equation does not 
has a non-trivial solution. 

• Algebraic general case. 

Consider the generic case 

Pi + P3, Pi -2P3, 3pi / 2p 3 . 
From the equations ( |6.28 - 3.27 ) the value of the constants are 



2pi 
2p3 - 3pi 



Pi ~ 2p 3 
2p3 - 3pi 



(6.33) 
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Now we derive by integration the ansatz functions 



f = a1n-o, g = b\n 



,2 ' 



(6.34) 



2« 



Thus the coframe is 



To 

The correspondent metric element is 



0°=(_) dt, &= — dx\ i = l,2,3 



ro 



2I> 



(6.35) 



ro 



la 



efcr= — dr- — (da^ + d^ + d^). 



ro 



46 



(6.36) 



Note that the parameters o and b are completely defined by Eq. (|6.33j ). 

Consider now how these proportional solutions related to the special solutions considered above. 

1) Pi = 

We obtain from the equations ( 6.33 ) 



a = 0, 



thus the solution (6.35) takes the form 



^ = dt, tf l = ( — ) dx\ i = 1,2,3. 



ds 2 = dt 2 -=( r -°Y(dr 2 + r 2 cm 2 



The correspondent metric is 



Using the following change of the radial coordinate 



we obtain the Minkowskian metric. Thus the Schwarzschild coframe yields a unique non-flat 
metric. 

2) pi = P3 

Substating this relain in fl6.33j ) we obtain 



6=1 



and the solution is 



ro 



dt, 0» 



ro 



dx l 



i = 1,2,3. 



This is a special case of ( 6.10| ). 
3) Pi + p 3 = 

The relation ( |6.33 ) gives 



5' 
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and the solution (6.35) is 



#> = (™)*dt, &=(?±)*dx i , i = 1,2,3. 



Thus the field equation has two non-trivial solutions. 

3) ps = 

The parameters of the cofrane are 

2 h 1 

a = — , o = — 

3 3 



and the solution is 



(-Ydt, &=(^)*dx i , i = l,2,3. 



Thus for P3 = the field equation has an unique solution. 

4) pi = -2p 3 

The parameters take the values 

1 h 1 

a = — , o = — 

2' 2 



and the solution is 



0° = ^cft, ^ = r -^dx\ i = 1,2,3. 



This is a special case of the conformal flat solution ( 6.31 ). 
5) pi = 2p 3 

The parameters a, b are non-defined and the proportional solution does not exist. 
The results of our analyse can be stated in the following theorem: 



Theorem 6.1: The field equation fiS.SJ has a non-trivial sphrically symmetric static solytion 
of the "diogonal" type for all values of the free parameters pi,ps exept of the case 

Pi = 2/5 3 . 

For 

Pi = P3 

exists a family of solutions, which depends on an arbitrary function of radial distance. 
For 

pi = -2p 3 

exists a family af conformal flat solutions, which parameterized by an arbitrary constant. 
For 

Pi + P3 = 

exist two non-trivial solutions. For all remain values of the parameters pi,P3 the solution is 
unique. 

The assymptotyc flat solution exist only for 

px = or pi + p 3 = 0. 



19 



A unique solution wich gives the Newtonian limit exists for 

Pi = 

and leads to the Schwarzschild metric. 



7 Metrics generated by proportional solutions. 



In the previous section it was shown that for generic values of free parameters of the theory the 
field equation ( |3.3| ) has a unique non-trivial static spherical-symmetric coframe solution of the 
form ( 6.35| ). This solution generates the metric ( 6.36j ). The physical meaning of these metrics 
require consideration of the topological properties of the manifold where the metrics can be 
realized. 

The camponents of the metric tensor ( 6. 36] ) have two coordinate singularities: one at the origin 
of coordinates r = and the second at the coordinate infinity r = oo. In order to clearify 
the nature of these singularities calculate the scalar curvature of the metric ( |6.36 ) Using the 
formula (6.1) we obtain 

4pi (a - b) 2 /r yb+2 



r= — , v - (7.i) 

pi + 2p 3 r 2 \r) 
Thus in the case 

26+1 >0 <==> H < -2 or H > - (7.2) 

P3 P3 3 

the scalar curvature in the origin r = is infinity. 
Accordingly in the case 

2b + 1 > <==> -2 < — < - (7.3) 

P3 3 

the scalar curvature at the second coordinate singularity r = oo is infinity. 
Note also the sign of the scalar curvature: 

R < <==> -2 < — < (7.4) 

P3 

and 

R > <==> f± > or P —<-1. (7.5) 

P3 P3 

Let us transform now the metric ( |6.36| ) to the Schwarzschild coordinates. In the spherical 
coordinates this metric takes the form 

ds 2 = (21 > ) 4 V_ (l-Y\dr 2 + r 2 dn 2 ). (7.6) 
\r / Vr / 
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In the case 

26+1/0 <==> 
we can use the transformation 



Pi 
P3 



7^-2 



f \ 2b 



, r \ 25+T 

> r ;•(, — 

,r , 



(7.7) 



(7- 



The metric ( |6.36 ) in the new (Schwarzschild) coordinates takes the form 

dr 2 
(26 + 1) 



dg2= , r_ytU dt 2 



r 2 dn 2 



(7: 



Note that that the hyper-spaces t = const is not flat, although 6 is constant. The four dimen- 
sional metric ( |7.9[ ) is spherically symmetric thus for analyzing the nature of the singularity it 
is enough to consider it's two-dimensional part "f — i" . The proper radial distance from an 
arbitrary point r\ to a singular point r* is 



I 



r* dr 




f* - h 


J fl 26+1 




26 + 1 



As for a proper time to attach the singular point r* by a radial null geodesic is 

df 



t 



/ /ro\2S+T 



|„m r m\ 
\ r * ~ r l I 



where 



26 + lJ Pl \rJ ' |26-2a + l|' 
26 - 2a + 1 2p 3 - 3pi 



(7.10) 



(7.11) 



m 



26 + 1 pi + 2p 3 

Thus we obtain the following regimes for the singularity in the origin of the coordinates: 



r = 


£i < -2 


-2< ^ < -i 


2 P3 3 


Pi v. 2 
P3 ^ 3 


New coordinates 


f = 


f = OO 


f = OO 


f = 


Scalar curvature 


R = oo 


= 


R = 


R = oo 


Proper distance 


finite 


oo 


oo 


finite 


Parameter m 


m < 


m > 


m < 


m > 


Proper time 


oo 


oo 


finite 


oo 



Analyses of this values gives the following regimes: 
^ < -2 and £i 



> ? 

P3 p3 3 

For these ranges of the parameters we obtain that in the point r 



the scalar curvature 

is singular. This point is located at a finite proper distance from an arbitrary point on 
the manifold but it takes an infinite time to reach this point by a radial null geodesic. 
Therefore for these values of the free parameters the origin of the coordinates is a physical 
singular point of the same type as the central point of the Schwarzschild metric. 
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_ 2 <£L< _| 

P3 -a 

The scalar curvature in the origin of coordinates is zero. This point is located on an 
infinite proper distance from an arbitrary point on the manifold and it takes an infinity 
time to reach this point by a null geodesic. Therefore in this case the point r = does 
not belong to any final portion of the space time. 

_I < £i < 2 

2 ^ p 3 ^ 3 

The origin r = is a point with zero scalar curvature. It is located in an infinity distance 
from every point but it can be reached by a final time. Thus the singular point is located 
outside of the null cone. 



As for the second coordinate singularity r = oo we have a following table of regimes: 



r = oo 


& < -2 

P3 


-2 < a < -\ 

P3 2 


2 Pi 3 


Pi v. 2 
P3 ^ 3 


New coordinates 


f = OO 


f = 


f = 


f = OO 


Scalar curvature 


R = 


R = oo 


i? = OO 


R = 


Proper distance 


oo 


finite 


finite 


oo 


Parameter m 


m < 


m > 


m < 


m > 


Proper time 


finite 


finite 


oo 


oo 



Thus we have the following regimes: 



£1 < -2 

P3 

The scalar singularity in the point r = oo is zero. This point is located on an infinity 
distance from every finite point on the manifold. Although it takes a finite time to reach 
this point from an arbitrary point on the manifold by a radial null geodesic. Thus the 
point r = oo is located outside of the null cone. 

_ 2 < £1 < -i 

P3 -a 

The scalar curvature is infinite thus the point r = oo is a physical singularity. This point 
is located on a finite distance from an arbitrary point and can be reached by a finite time. 
The singularity is similar to a vertex of the cone - conical singularity. 

_I ei. <r ^ 

2 ^ p 3 ^ 3 

For these parameters the point r = oo is, also, a physical singularity. It can be however 
be reached by a infinite time nevertheless it is located on a finite proper distance. Thus 
this singularity is similar to the central point of the Schwarzschild metric. 

Pi ^ 2 
P3 ^ 3 

The proper distance and the proper time is infinite. Therefore the point is located outside 
of any finite area of space-time. 
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8 Conclusions and discussion. 



We study the general quadratic Lagrangian of pure teleparallel theory. This Lagrangian is a 
linear combination of three covariant and global SO(l, 3) invariant terms with free dimensionless 
coefficients pi,p2,P3- The field equations of the theory is studied by a "diagonal" ansatz which 
is a subclass of a general spherical-symmetric Einstein-Mayer ansatz. The restriction is taken 
for simplification the calculation. Thus we study a Lagrangian that depends only on two free 
parameters p\,pz- 

We obtain a formula for scalar curvature of a pseudo-Riemannian manifold with a metric 
constructed from the solution of the field equation. This formula shows that the sign of the 
scalar curvature depends only on a relation between the parameters p\ and ps- The scalar 
curvature vanishes only in a case r± = which corresponds to the teleparallel equivalent of GR. 
We obtain exact solutions for all possible values of free parameters pi,pz- It is shown that the 
unique solution with a Newtonian limit is the Schwarzschild solution. Thus the Yang- Mills-type 
term of the general quadratic Lagrangian should be rejected. 

Note that all the results valid only in the case of a restricted "diagonal" ansatz. Thus it is 
important to obtain solutions for following questions: 

• Which physical sense can be given to the various exact solutions? 

• How the results depends on the non-diagonal terms? 

• How the scalar curvature inequalities can be generalized to the case of three free param- 
eters Pl,p2,P3? 
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A Calculations with the diagonal ansatz 



We study the field equation (3.3) with the diagonal coframe ansatz 

d° = e f dt, d l = e 9 dx l , i = l,2,3. 

We will use the following notations 

(V/ ■ Vg) = + /2S2 + f 3 g 3 
(V/) 2 = A 2 + /| + ft 

A/ = fll + / 2 2 + fss 
Calculating the first derivatives of the diagonal coframe we obtain 



(A.l) 
(A.2) 
(A.3) 





*dtf° = -e- 9 (/n9 23 - / 2 tf 13 + / 3 tf 12 ) 


d-d 1 = -e- 9 (g 2 D 12 + g 3 d 13 ) 




d-d 2 = e~ 9 ( gi -d 12 -g 3 -d 23 ) 


*dtf 2 = - e - s (pitf 03 -p 3 i? 01 ) 


di9 3 = e' 9 (gid ri + 321? 23 ) 





A.l Temporal components 

The first leading term of (^^) takes the form 

(1) L = d*dtfo = -e" 29 (A/ + (V/-Vg))# 123 , 

Thus 



(i) 



Lo = -e- 29 (A/ + (V/-V ff ))* 1 9° 



(A.4) 
(A.5) 



As for the second leading term of (3.3) 

(3) L = A d * (tfo A dtf 6 ) = -i? 1 A d * A d-d 1 ) - -d 2 A d * A di9 2 ) — # 3 A d * A dtf 3 ) 
= A d * (-e- 9 )( ff2 tf 012 + g 3 d 013 ) -d 2 Ad* {e- 9 )( gi i) m2 - g 3 d 023 ) 

—d 3 Ad* (e- 9 )( ff ij? 013 + g 2 d 023 ) 
= d 1 A d(e" 9 )(<72J? 3 - sstf 2 ) — i? 2 A d(e~ 9 )( 5 n9 3 - ^tf 1 ) + # 3 A dfe -9 )^!? 2 - gid 1 ) 
= & 1 A d(g 2 dz - g 3 dy) - d 2 A d(gidz - g 3 dx) + -d 3 A dfoidy - g 2 da;) = 2e~ 29 Ag-d 123 

Thus 



< 3 >L = 2 e - 29 A g *tf° 



(A.6) 
(A.7) 



The first quadratic term is 

(1) Qo = eoJ {d-d b A *d-d b ) = e J (d#° A *d-d° - dtf 1 A *dtf 1 - d-d 2 A *di9 2 - d-d 3 A *dtf 3 ) 

= e~ 29 e J ((fid 01 + f 2 d 02 + f 3 d 03 ) A (fid 23 - f 2 d 13 + f 3 d 12 ) - (g 2 d 12 + g 3 d 13 )(-g2d 03 + g 3 d 02 ) + 
(gid 12 + g 3 d 23 )(~gid 03 - g 3 d 01 ) - ( 9l d 13 + 9 2d 23 ){gid 02 - g 2 d 01 ) 



e - 29 {V 2 f + 2V 2 g)$ 



(A.8) 
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Thus 



(1) Qo = e- 29 ((V/) 2 +2(V<;) 2 )* 1 ? ,) 



The second quadratic term is 

(2) Qo = (eoJdtf 6 ) A *<M b = (e Jdtf°) A *dD 

= e-^ihi) 1 + / 2 tf 2 + / 3 tf 3 ) A (/itf 23 - / 2 tf 13 + / 3 tf 12 ) = e- 29 V 2 /^ 123 

Thus 



(2) 



Qo = e- 29 (V/) 2 *^ 



The next non-zero quadratic term is 



(6), 



»o = dti b A *{& A dd b ) = -d-d 1 A *(#° A dtf 1 ) - dtf 2 A A di? 2 ) - di? 3 A A dtf 3 ) 



-2s 



-29 



((<? 2 tf 12 + fl3 13 ) A *(S2tf° 12 + 53 tf 013 ) + ( gi {> 12 - 53 tf 23 ) A *( 5l tf 012 - S3 

+ ( 5 i^ 13 + 32^ 23 ) A *( Sl tf 013 + ff2 tf 023 )) 
((52i9 12 + Satf 13 ) A (<? 2 tf 3 - g 3 $ 2 ) + ( gi $ 12 - <? 3 tf 23 ) A ( S itf 3 - grf 1 ) 

+( Sl tf 13 + 52 tf 23 ) A (- 5l tf 2 + gjtf 1 )) = -2e- 29 V 2 51 9 123 



Thus 



(6) 



Qo = -2e-' 3 (Vgy* 



As for the next non-zero quadratic term ^Qo we have 
(7) Qo = e J (tf c A dtf b A *{dd c Ad b f) 

= 2e J (-/A d-d 1 A A#)-t)°A d# 2 A *{dd° A tf 2 ) - A dtf 3 A A tf 3 ) 

A dtf 2 A *(dtf 1 A j? 2 ) + 1? 1 A dtf 3 A A # 3 ) + $ 2 A dtf 3 A *(dtf 2 A #)) 

= 2e- 2 «e J (( 52 tf 012 + 53 #' 13 ) A (/ 2 tf 3 - / 3 tf 2 ) + ( 5l tf 012 - gtf 023 ) A (/^ 3 - /stf 1 ) 
+( gi # 013 + ff2 tf 023 ) A (-f^ 2 + M 1 ) - g^ 12i A g^° - <? 2 tf 123 A S2 tf° - ffl tf 123 A 
= 2e~ 2s (2V/ ■ Vg + V 2 <?)# 123 



Thus 



(7) Q = 2e- 29 (2(V/-V ff ) + (V ff ) ; 



The last non-zero quadratic term ^Qo is 

(8) Qo = (eoJdtf 6 ) AtfcA *(dtf c A S b ) 

= -(e Jdtf°) A (V A *{d^ A 0°) + # 2 A *(dtf 2 A + t? 3 A *(d# 3 A tf )) 

= -e" 2 ^/^ 1 + / 2 tf 2 + /a,? 3 ) A (V A *(<? 2 #' 12 + 53 #' 13 ) — t? 2 A *(<7i0° 12 - AS* 023 ) - 

tf 3 A*( Sl tf 013 + S2 O) 
= -e-^/itf 1 + / 2 tf 2 + / 3 tf 3 ) A ( 52 t? 13 - 53 tf 12 - <?itf 23 + 53 tf 21 + gi t? 32 - g 2 tf 31 ) 



2e- 29 (V/- Vff)^ 1 



Thus 



(8), 



) = 2e- 29 (V/-Vff)*^° 
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A. 2 Spatial components 



The first leading term takes the form 

(1) Li = d * d&i = d(e f (-g 2 dt Adz + g 3 dt A dyjj 

= -e- 2a (( Sl2 + /i<?2)#' 13 + (322 + 033 + / 25 2 + fsgsW 023 - (013 + fig 3 )ti 012 ) 
= e- 2g ((A3 + V/ • Vg - 311 - figi) * ti 1 - (312 + /132) * tf 2 - (313 + figs) * V s ) 

Thus 



(i) 



Li 



,-29 



A 9 + (V/ • V3) W + ftfc + * 



The second leading term is 

(3) L 1 =D b Ad*(-d 1 Add b ) 

= Ad* (d 1 A d$°) + tf 2 A d * (d 1 A d-d 2 ) + i9 3 A d * (tf 1 A d-d 3 ) 

= -#° A d * (e" 9 )^ 012 + / 3 ^° 13 ) - t? 2 A d * (e- 9 ) 53 i? 123 ) + d 3 A d * (e^ 9 ^ 123 ) 

= -d° A d{e- 9 )(f 2 d 3 - f 3 d 2 ) - d 2 A d(e- 9 g 3 d°) + d 3 A d{e~ B g 2 'd ) 

= -d° A d(f 2 dz - f 3 dy) — -d 2 A d(e f ~ s g 3 dt) + d 3 A d{e f - 9 g 2 dt) 

= -d° A ((/ 2 2 + f 33 )dy A dz + f 12 dx Adz - f 13 dx A dyj 

-d 2 A e f - 3 ([313 + (/1 - 31)33]^ A dt + [333 + (/s - gz)gs]dz A dt) 

+d 3 A e f - 9 ([g 12 + (/1 - gi)g 2 ]dx A dt + [322 + (/ 2 - 32)32]*% A dtj 

= e~ 29 ((-A/ - A3 + (V3) 2 - V/ • V3 + /n + 311 - 3i + Affi) * V 1 

+ (/l2 + 312 + /l32 - 3132) * $ 2 + (/l3 + 313 + /l33 - 3153) * ^ 3 ) 

Thus 



(3) i, = e- 29 (^Af + A3 - (V3) 2 + (V/ • Vg)) Vik + f ik + g ik - 9l g k + f t9l ?J * ti k 



Using the calculations of ^Qo we obtain the first quadratic term as 



= eiJ (d$ b A *d$ b ) = e- 29 ((V/) 2 + 2(V3) 2 )^ fe * # k 



(A.18) 



(A.19) 



(A.20) 



(A.21) 



(A.22) 



As for the second quadratic term 

(2) Qi = (ejdtf 6 ) A *dti b = (eiJdtf ) A *dD° - (eddtf 1 ) A *d-d 1 - (eiJdd 2 ) A *d-d 2 - (ejdtf 3 ) A *<M 3 
= e- 29 (-/A (/h? 23 - f 2 d 13 + f 3 d 12 ) - ( 52 tf 2 + 33^ 3 ) A (-32^ 03 + 33t?° 2 ) 

-3ii9 2 A (3it? 03 - g 3 -d in ) - gi-d 3 A ( fl n?02 - 3 2 tf 01 )) 
= e" 29 ( - (/1 + 32 2 + gt + 2gl)d 023 + (/i/ 2 + 3i3 2 )tf 013 - (/1/3 + 3i33)^" 12 ) (A.23) 



Thus 



(2) 



Q* = e- 29 {{Vgfr)ik - {fif k + g l9k )) * d k 



(A.24) 
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The next non-zero quadratic term is 

(6) Qi = d-d b A *(#i A d-d b ) = ~d-d° A *(-d 1 A d-d ) + dd' 2 A *(-d 1 A d-d 2 ) + d-d 3 A ♦ (tf 1 A ctt> 3 ) 
= e- 29 ((/i^ 01 + /a/" + / 3 t? 03 ) A *(/ 2 tf° 12 + jatf 013 ) - 

( 5 itf 12 - 53 ^ 23 ) A *g :i i) 123 + (ffi^ 13 + 32# 23 ) A * 521 9 123 ) 
= e- 29 ((/xtf 01 + / 2 tf 02 + / 3 t? 03 ) A (/ 2 tf 3 - / 3 tf 2 ) - (gid 12 - g z d 23 ) A S3 ^° + (gi$ 13 + 52i? 23 ) A * S2 tf°) 



e- 29 ((/I + /? + si + Sl)^ 023 + (/1/2 + Si<? 2 )tf° 13 - (/1/3 + Siffs)^ 012 ) 



Thus 



(6)/ 



>i = -e 



-2g 



((V/) 2 + ( V 5 ) 2 ) Tfcfc +fif k + grg k 



(A.25) 
(A.26) 



Using the calculations above for ^Qo we obtain for the next non-zero quadratic term 

(7) <2i = eiJ (# c A d$ b A *(dr A = -2e-' 29 (2Vf ■ V5 + V 2 g) * tf 1 

Thus 



(7) Q, = 2e- 29 (2(V/ • V5) + (\7g) 2 )r, lk * 



(A.27) 



As for the last non-zero quadratic term 
(8) Qo = (ejdtf 6 ) A tf c A *(dt? c A 

(V A A + tf 2 A *(<M 2 A + i? 3 A *(<M 3 A tf )) 

(-/A A tf 1 ) + tf 2 A *(d'i? 2 A i? 1 ) + # 3 A *(dtf 3 A i? 1 )) 

(-/A A tf 2 ) + i? 1 A ♦ (dtf 1 A i? 2 ) + $ 3 A *(dtf 3 A tf 2 )) 

(-/A A tf 3 ) + i? 1 A ♦ (dtf 1 A i? 3 ) + i? 2 A *(<M 2 A tf 3 )) 

(V A *( 52t ? 012 + g 3 $ 013 ) — -d 2 A *(g^ 012 - g 3 # 023 ) - tf 3 A *( 9l $ 013 + 32^° 23 )) 



= -(ejdtf ) A 
+ (eiJ(M 1 ) A 
+ (eiJ<M 2 ) A 
+(eiJdtf 3 ) A 



,-2g 



+ (521? 2 + fl 3 tf 3 ) A /A *(/2i? 012 + / 3 ^° 13 ) + # 2 A 



-29 



3 + J? 3 A *g 2 tf 123 ) 
+ 5 i -d 2 A (&° A *(/itf 012 - / 3 tf 023 ) + i? 1 A * 33 1? 123 — i9 3 A *<? 3 tf 123 ) 

+51 # 3 A (V A *(/itf 013 + / 2 tf 023 ) - ^ A * 52 tf 123 + d 2 A *ffitf 123 )) 
/11? A (V A (52^ - g 3 ti 2 ) - -d 2 A ( 5 itf 3 - 531? 1 ) + -d 3 A ( 3 ii9 2 - 521? 1 )) 

+ (<?2# 2 + gs-d 3 ) A A (/ 2 t? 3 - /ai? 2 ) + -d 2 A 5 3#° + # 3 A ff2 l?°) 
+ 5 itf 2 A A (/itf 3 - ft-d 1 ) + -d 1 A gs-d" - V 3 A fl3 0°) 

+5it? 3 A (d° A (-/itf 2 + / 2 tf x ) — t? 1 A 52 i9° + i9 2 A 51 



-29 



(4/151 + / 25 2 + / 3 53 + 2g! + gl + <?3 2 )tf 023 + (2/152 + 51/2 + 5i<? 2 )tf 013 



-(2/153 +51/3 +5i5 3 )^ 012 ) 



Thus 



(8) 



Qi = e 



-2s 



(((V/ • Vg) + (V5) 2 )r; lfc - 2f i9k - 5,/ fc - 5,5^ * -d k 



(A.28) 
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